Axisymmetric solitary waves on the surface of a ferrofiuid 
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We report the first observation of axisymmetric solitary waves on the surface of a cylindrical 
magnetic fluid layer surrounding a current-carrying metallic tube. According to the ratio between 
the magnetic and capillary forces, both elevation and depression solitary waves are observed with 
profiles in good agreement with theoretical predictions based on the magnetic analogue of the 
Korteweg-deVries equation. We also report the first measurements of the velocity and the dispersion 
relation of axisymmetric linear waves propagating on the cylindrical ferrofiuid layer that are found 
in good agreement with theoretical predictions. 

PACS numbers: 47.35.-i,47.65.Cb,47.35.Fg 



Solitary waves or solitons are localized nonlinear waves 
that propagate almost without deformation due to the 
balance between the nonlinearity and the dispersion. 
Since the first observation of a solitary wave on the free- 
surface of water by Russel [l|, and its interpretation us- 
ing the Korteweg-de Vries equation (KdV) [2], it has 
been shown that the KdV equation describes a large 
class of solitons observed in various situations: acous- 
tic waves on a crystal lattice, plasma waves, hydrody- 
namics internal or surface waves, elastic surface waves, 
and waves in optical fibers Q. Most of them involve 
a localized elevation disturbance propagating within a 
quasi-one-dimensional plane system. Observations of ax- 
isymmetric solitary waves governed by the KdV equation 
are scarce and mainly concern waves in rotating flu- 
ids confined in a tube or on vortex lines. More recently, 
Bashtovoi et al. derived a KdV equation with an axisym- 
metric solitary waves solution propagating on the surface 
of a cylindrical magnetic fluid layer submitted to a mag- 
netic field 0, Without gravity, the stability of the 
cylindrical magnetic fluid layer is governed by the ratio 
between the magnetic force and the capillary one. Ac- 
cording to its ratio, axisymmetric elevation (hump-like) 
or depression (hole- like) solitary waves are predicted with 
a subsonic or supersonic velocity 0, @ . To our knowl- 
edge, direct observation of axisymmetric magnetic soli- 
tary waves have never been reported. 

In this Letter, we report the first observation of ax- 
isymmetric solitary waves on the surface of a cylindri- 
cal ferrofiuid layer submitted to an azimuthal magnetic 
field. Depending on the strength of the field, elevation or 
depression solitary waves are observed on the ferrofiuid 
surface. A ferrofiuid is a stable suspension of nanometric 
magnetic particles diluted in a carrier liquid (water or oil) 
that responds to an external applied magnetic field 0, [1| . 
Although the solitary waves are damped by viscous dis- 
sipation, we have shown that they keep the self-similar 
profile given by the solution of the KdV equation on a 
propagation length larger than their typical scale. More- 
over, we also report the first measurement of the velocity 
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FIG. 1: (Color online) Experimental setup. 



and dispersion relation of axisymmetric magnetic linear 
waves in this system in good agreement with the theo- 
retical predictions 0. 

The experimental setup is shown in Fig. [TJ It consists 
of a cylindrical copper tube (50 cm in length, a = 1.5 
mm in outer radius and 0.5 mm in thickness) placed in 
the middle of a rectangular plexiglas container (40 x 40 
mm 2 side and 30 cm length) crossing both container end- 
sides in the center through hollow waterproof screws. A 
dc electrical current, /, in the range 0-100 A is applied 
to the cylindrical conductor by means of a power supply. 
The current generates circular magnetic field lines around 
the tube with a radial decreasing amplitude. The corre- 
sponding radial magnetic force stabilizes a ferrofiuid layer 
of outer radius 2.9 < R < 5 mm, around the tube. For 
the results reported below, R — 3.8 mm for linear waves 
and R = 3.3 mm for solitary waves. The ferrofiuid used 
is an ionic aqueous suspension synthesized with 12.4% 
by volume of maghemite particles (Fe203 ; 7 ±0.3 nm in 
diameter) [l(| • The properties of this magnetic fluid are: 
density, p = 1534 ±1 kg/m 3 , initial magnetic susceptibil- 
ity, Xi = 0.75, magnetic saturation M sat = 36 x 10 3 A/m, 
and estimated dynamic viscosity 1.4 xl0~ 3 Ns/m 2 . To 
avoid gravitational effects, the whole container is filled 
with Freon (C2CI3F3), a nonmiscible transparent fluid 
with a density, p e — 1581 ± 1 kg/m 3 , close to the fer- 
rofiuid one. The surface tension between the ferrofiuid 
and freon is 7 = 5.5 x 10~ 3 N/m. A water cooling inside 
the tube drains off the Joule dissipation of the current- 
carrying electrodes (their contact resistance being 13.8 
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mf2, that is a dissipated power of roughly 140 W for 
/ = 100 A). The metallic tube is tightened to avoid par- 
asitic vibration. To wit, one end of the tube is threaded 
to fix it with a nut, and a small chuck is used at the 
other end to tighten it. Surface waves are generated on 
the fcrrofluid surface by the horizontal motion of a con- 
centrical plexiglas tube, 9 mm (resp. 3 mm) in outer 
(resp. inner) diameter driven by an electromagnetic vi- 
bration exciter. The wavemaker is driven sinusoidally (in 
a frequency range from 0.5 to 10 Hz with a maximal am- 
plitude of 2 mm) to study linear waves, or impulsively 
(typical duration of 0.05 s) to study solitary waves. Note 
that the wavemaker end is made of copper in order to in- 
crease the ferrofluid wetting. Axisymmetric waves prop- 
agating on the cylindrical ferrofluid layer are visualized 
with a high-resolution camera (Pixclink 2208 x 3000 pix- 
els) located above the container, and are recorded with a 
25 Hz (resp. 44 Hz) sampling for linear (resp. solitary) 
waves. 

The magnetic induction generated by the carrying- 
current tube is up to 30 G at 100 A at a distance r = 8 
mm from the tube axis, that is the z-axis of the (r, 9, 
z) cylindrical coordinate system. The magnetic induc- 
tion being orthoradial B—[B r = 0, Bg = B, B z = 0], i.e. 
throughout tangential to the free surface, the Rosensweig 
magnetic surface instability is absent 0, [ll| . B is mea- 
sured with a transverse Hall probe via a gaussmeter (Bell 
5100) as a function of the current and the distance r from 
the tube in agreement with the usual law B = p^I / '(2wr) 
for r > a where po = 47r x 10 H/m is the magnetic 
permeability of the vacuum. The corresponding mag- 
netic body force, F mag = — p xl 2 / {4:Tr 2 r 3 )e r , is radial 
towards the z-axis. Since gravity is negligible (p e ~ p), 
this magnetic force stabilizes a uniform axisymmetric 
layer of a magnetic fluid of outer radius R as soon as 
the capillary force by volume F cap ~ ijR 2 is small 
enough. Both magnetic and capillary effects are then 
compared by the dimensionless magnetic Bond number 
Bo m ee F mag /F cap = p xI 2 /(4:ir 2 jR). 

Assuming no gravity and a thin tube radius (a <C R), 
the dispersion relation of inviscid axisymmetric linear 
waves propagating on a magnetic fluid surface reads 
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where uj = 2?tf is the angular frequency and k ee 2tt/X 
the wavenumber, /„ and K n being, respectively, the mod- 
ified Bessel functions of first and second kind of order n 
(their ratio being a positive increasing function of k). 
When Bo m < 1, the capillary effects are greater than the 
magnetic ones, and an instability occurs [uj 2 < in Eq. 
(fl|)]: the cylindrical ferrofluid layer is unstable to dis- 
turbances whose wavelengths A > 2-7ri?/-\/Bo m — 1, and 
breaks up into a string of connected drops 0, [lj| . This 
is the magnetic analogue of the surface-tension-driven 




FIG. 2: (Color online) Dimensionless dispersion relation of 
linear cylindrical waves for various applied currents / from 40 
to 100 A corresponding to Bo m = 1.85 (x), 2.69 (+), 4.17 (<), 
4.89(A), 5.67 (0), 6.51 (o), 7.4 (v), 8.36 (0), 9.37 (*), 10.44 
(□), and 11.57 (>). Solid line corresponds to the theoretical 
prediction of Eq. {]]). Top inset: Frequency / as a function 
of the wavelength A for different / with a 10 A step. Dashed 
lines are from Eq. {]]). Bottom inset: Snapshot of linear waves 
(/ = 3.5 Hz, Bo m = 6.5 - wavemaker is visible on the left-hand 
side - 10 cm size window). 



Rayleigh-Plateau instability when a thin cylindrical jet 
of a usual fluid breaks into a set of drops [l3j]. When 
Bo m > 1, one has uj 2 > in Eq. ([I]): the cylindrical 
layer of ferrofluid is stable whatever the wavelength dis- 
turbance, and axisymmetric linear waves can propagate 
on its surface. 

We first measure the dispersion relation of such linear 
waves. The wavemaker is driven sinusoidally in order to 
generate surface waves at the interface between the freon 
and the fcrrofluid. A typical snapshot of such axisymmet- 
ric linear waves is shown in the bottom inset of Fig. [2] 
The top inset of Fig. [5] shows the wavelength, A, of surface 
waves in response to the forcing frequency, /, for differ- 
ent applied currents, /, that is for different Bo m ~ I 2 
ranging from 1 to 12. A is found to decrease with in- 
creasing frequency whatever Bo m . When expressed in 
the rescaled variables uj 2 / [(Bo m - 1 + (kR) 2 ) ~//(pR 3 )] 
and kR, all these data collapse on one single master curve 
(solid line) predicted by Eq. ([T]) (see Fig.rj]). Note that no 
adjustable parameter is used when comparing the data 
and the theoretical dispersion relation of axisymmetric 
magnetic surface waves. 

Using the expansion of the modified Bcsscl functions 
— fg lM > the dispersion relation of Eq. (JXJ) in 
a long- wavelength limit (kR <C 1) reads @ 
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Co being the velocity of linear waves for Bo m > 1 . Note 
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FIG. 3: (Color online) Velocity of linear axisymmetric mag- 
netic waves as a function of rescaled Bond number. Experi- 
mental (♦) and theoretical [solid line of Eq. linear wave 
velocity cq. Dashed line: theoretical solitary wave velocity c 
of Eq. (JS| with Ao = 0.5 mm, R — 3.8 mm (see text). Dash- 
dotted lines correspond to Bo m = 3/2 and 9. Inset: Rescaled 
dispersion relation uj /(co/R) vs. kR. Same symbols as Fig. 
[2] Solid line has a slope 1. 



that both Co and the sign of the dispersive term ~ k 3 in 
Eq.© depend on Bo m . In order to extract the velocity 
of linear waves, the inset of Fig. [3] displays the previous 
data in variables ui/(cq/R) and kR. As expected, for 
small kR, all the data collapse on a single linear curve 
of slope 1. Note that for larger kR and for Bo m ~ 1, a 
departure from the prediction of the linear term of Eq. 
@ is observed since the dispersive effects become im- 
portant [i.e., the fc 2 /(Bo m — 1) term in Eq. ©]. For 
kR <C 1, the slope of each curve w vs. k thus gives, for 
each Bo m value, a direct measurement of the velocity of 
linear waves. These values are plotted in Fig. [3] and are 
found in rough agreement with the theoretical velocity of 
Eq. © with no adjustable parameter. To our knowledge 
this is the first measurement of the velocity and disper- 
sion relation of axisymmetric magnetic linear waves on 
the surface of a ferrofluid. 

Let us now focus on axisymmetric magnetic solitary 
waves. First, let us assume no viscosity and no gravity. 
In the long- wavelength limit (kR <C 1), the dispersion is 
small and the linear wave velocity is Co- When the inter- 
face deflection A(z,t) is also small, such that nonlinear 
effects have the same order of magnitude as dispersive 
ones, it is governed at the leading order by a magnetic 
analogue of the Korteweg-de Vries equation 0] 
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the nonlinear coefficient, 



j w — — the dispersive coefficient, and 



Bo m > 1. The axisymmetric magnetic solitary wave so- 
lution of Eq. (3) reads 0, 
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and L is the length scale of the solitary wave. Equations 
(HJ and ([5]) show that there exists a continuous family 
of soliton solutions with parameter Ao (the extremum 
amplitude of the wave). Since a, (3, c and L depend on 
Bo m , an elevation (A > 0) or depression (Ao < 0) soli- 
tary wave is predicted that propagates with a supersonic 
(c > c ) or subsonic (c < Co) speed. All the possible 
solutions are summarized in Tab. HI 
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TABLE I: Properties of axisymmetric solitary wave solution 
of Eqs. and ([5]) according to the magnetic Bond number. 



We have performed a study of axisymmetric solitary 
waves on the surface of the ferrofluid layer (R — 3.3 mm) 
around a copper tube carrying current in the range 60 to 
110 A (4 < Bo m < 14). We impulsively drive the shaker 
to generate solitary waves: the wavemaker is pushed for- 
ward to generate a pulse on the fluid interface leading 
to either an elevation or a depression pulse according 
to the value of Bo m . The interface deflection A(z,t) is 
detected from the images recorded by the camera using 
standard ImageJ binarization and edge detection pro- 
cesses. The profile is displayed in Fig. @Ji for a depres- 
sion pulse (Bo m = 8.3) and in Fig. |4p for an elevation 
pulse (Bo m = 10.5). Both recordings are in good agree- 
ment with the profiles of elevation and depression KdV 
solitary waves given by Eq. (QJ. Note that once Ao is 
known the theoretical profile as well as the velocity of 
the solitary wave given by Eqs. ((4]) and ([5]) do not in- 
volve any adjustable parameter. Those isolated pulses 
involve typical amplitude and size that are in the range 
of validity required for the derivation of Eq. ©, that is, 
corresponding to small dispersion (L 2 3> R 2 ), and small 
nonlincarities (\Aq\ <C R), both of same order of mag- 
nitude (R 3 ~ |j4o|£ 2 )- Note that no solitary wave has 
been observed for Bo m < 4 since its predicted amplitude 
(Ao ~ few mm for L ~ 1 cm) is too large compared 
to our tube radius R to have small nonlinearities. In- 
set of Fig. QJi (resp. Fig. 0J)) shows the profile of the 
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FIG. 4: (Color online) (a) Axisymmetric profile of a de- 
pression magnetic solitary wave for Bo m = 8.3 (I = 80 
A), t = 0.33 s, (b) elevation magnetic solitary wave for 
Bo m = 10.5 (I = 90 A), t = 0.32 s, centered on its extremum. 
Solid lines are the theoretical profiles of KdV solitons derived 
from Eq. @ with no adjustable parameter. Insets: Profiles 
of solitary waves at different times [(a): 0.16, 0.33, 0.47, 0.62, 
and 0.74 s; (b): 0.17, 0.32, 0.47, 0.58, and 0.65 s] during its 
propagation over 10 times its typical size L ~ 4 mm. Initial 
amplitudes Aq — ±0.5 mm. R — 3.3 mm. z-axis origin is 
located on the wavemaker. 



depression (resp. elevation) pulse recorded at different 
times corresponding to a total propagation distance up 
to 10 times its typical size. The recorded profiles are in 
good agreement with the KdV magnetic solitary wave 
all along the propagation. Note, however, that for far- 
ther distances the cumulative effect of dissipation leads 
to small amplitudes that are hardly measurable by the 
camera (0.04 mm/pixel). For both the elevation and de- 
pression solitary waves, dissipation leads to an extremum 
amplitude Aq(z) that decreases linearly with the propa- 
gation distance z. By rescaling all the profiles displayed 
in each inset of Fig. 0] with the variables |A(z)/Ao(z)| 
vs. z/L, all the data lie on a single curve predicted by 
Eq. ((Ij) (not shown here). This means that the pulse 
keeps a self-similar shape over a distance up to 10 times 
its typical size and is in a good agreement with the pro- 
file derived from KdV magnetic equation. Finally, the 



solitary wave velocity, c, is measured all along its propa- 
gation by the successive locations of amplitude extrema, 
|j4 (z)|, at different times. We find that c ~ cq (~ few 
cm/s) with a dependence on Bo m roughly comparable to 
that predicted by Eq. (T2]). Note that c is predicted by 
Eq. (|5|) to slightly depend on Aq with a correction with 
respect to c up to 10% when Bo m ~ 10 (see dotted line 
on Fig. [3]). Since our velocity measurement accuracy is 
6%, we cannot thus discriminate from a subsonic to a 
supersonic solitary wave as predicted in Table I (see also 
Fig. [3D . 

In conclusion, we have reported the first observation 
of depression and elevation axisymmetric solitary waves 
on the surface of a cylindrical magnetic fluid layer and 
found that their shapes are in good agreement with the 
ones predicted from the axisymmetric KdV solitary wave 
solutions. A possible extension of this work would be the 
study of the collisions between these new solitary waves. 
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